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Abstract 

We explore the synchronization behavior in interdependent systems, where the one-dimensional 
(ID) network (the intranetwork coupling strength Ji) is ferromagnetically intercoupled (the 
strength J) to the Watts-Strogatz (WS) small-world network (the intranetwork coupling strength 
Jn). In the absence of the internetwork coupling (J = 0), the former network is well known not 
to exhibit the synchronized phase at any finite coupling strength, whereas the latter displays the 
mean-field transition. Through an analytic approach based on the mean-field approximation, it is 
found that for the weakly coupled ID network (Jj <C 1) the increase of J suppresses synchrony, 
because the nonsynchronized ID network becomes a heavier burden for the synchronization process 
of the WS network. As the coupling in the ID network becomes stronger, it is revealed by the 
renormalization group (RG) argument that the synchronization is enhanced as J\ is increased, im- 
plying that the more enhanced partial synchronization in the ID network makes the burden lighter. 
Extensive numerical simulations confirm these expected behaviors, while exhibiting a reentrant be- 
havior in the intermediate range of J\. The nonmonotonic change of the critical value of Jn is also 
compared with the result from the numerical RG calculation. 
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Synchronization phenomenon is the one of the most fascinating collective 
emergent behaviors abundantly found in natural and artificial systems. The on- 
set of synchronization occurs when the differences of individual oscillators are 
overcome by the strong coupling among elements. We investigate the coupled 
system of two networks, one is synchronizable at finite coupling strength while 
the other is not, aiming to answer the question of what happens as the inter- and 
intra-couplings are varied. For the weak intracouplings of the nonsynchronizable 
network, both our analytic and numerical results show that the stronger internet- 
work coupling hinders the synchronization in the synchronizable network since 
the nonsynchronized oscillators in the other network work as heavier burdens 
for the oscillators in the synchronizable network to carry to synchronize. On 
the other hand, as the intracoupling strength in the nonsynchronizable network 
becomes larger, partially synchronized groups of more oscillators are formed, 
which in turn help the oscillators in the synchronizable network to become uni- 
fied as synchronized clusters. In the intermediate regime where the intra- and 
inter-network couplings are in the same order of magnitude, numerical results 
show a reentrant behavior in the synchronization phase diagram. 



I. INTRODUCTION 



Synchronization as a collectively emergent phenomenon in complex systems has attracted 
much interest thanks to the abundance of examples in nature l|-|4j]- In existing studies, it 
has been revealed that the topology of interaction plays an important role in synchroniz- 
ability . In particular, it is now well known that coupled phase oscillators described 
by the celebrated Kuramoto model [4j exhibit various universality classes depending on di- 
mensionality and the topological structure of networks js il]. It has also been found that 
the lower critical dimension for the frequency synchronization in the regular <i-dimensional 
Kuramoto model is d = 2, while the corresponding lower critical dimension for the phase 



synchronization transition with the spontaneous 0(2) symmetry breaking is d — 4 [10|, 111 ]. 
Moreover, numerical studies H[ have implied that systems of d > 4 belong to the mean-field 
(MF) universality class as shown in globally coupled oscillators. On the other hand, for the 
Kuramoto model in complex networks, various results have been reported: For random and 
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small- world networks [jjl H], it has been observed that the MF transition exists js, 7], and 
for scale-free (SF) network with a power-law degree distribution p ~ A;" 7 , where k stands 
for degree, 7-dependent critical exponents have been found via a MF approximation and 
numerical investigations {9]. 

In the present work, we couple the one- dimensional (ID) network and the Watts-Strogatz 
(WS) small- world network [13| (see Fig. [T]), and investigate the synchronizability of Ku- 
ramoto oscillators in the composite two coupled networks. Our research focus is put on the 
effect of internetwork coupling between the two networks belonging to different universali- 
ties, i.e., the MF universality for the WS network and the absence of the synchronous phase 
for the ID regular network. The physics of coupled interdependent networks are not only 
interesting in the pure theoretical point of view, but it also can have practical applicabil- 
ity since these interdependent network structures can be found ubiquitously. For example, 
electric power distribution in the power grid is strongly interwoven with the communication 
throu gh I nternet, and thus spreading of failures in one network affects failures in other net- 



work 
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15]. This coupled system has been studied within the framework of percolation 



with the strength of internetwork coupling varied; strong coupling between networks yields a 
first-order transition, while in weakly coupled networks a giant cluster continuously vanishes 
at the critical point [if]]. In Ref. 17, the epidemic spreading behavior has been studied in 



the coupled networks of the infection layer and the prevention layer. We believe that the 
study of the collective synchronization in interdependent networks can also be an important 
realistic problem in a broader context: Imagine that each agent in a social system has two 
different types of dynamic variables and that the interaction of the one type of variable 
has different interaction topology than the other variable. We emphasize that our study 
of the synchronization in coupled networks is worthwhile because it could be applicable to 
investigate social collective behaviors in interdependent networks. Another interesting ex- 
ample of the interdependent system can be found in the neural network in the brain: The 
cortical region is coupled with the thalamus [l8|, and the thalamocortical interactions might 
be interpreted as the intercoupling between the cortical area and thalamus. 
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FIG. 1: (Color online) Two coupled interdependent networks: The one-dimensional regular network 
(the upper one) with the periodic boundary condition is coupled to the Watts-Strogatz small-world 
network (the lower one). Both networks have the average degrees four. 

II. MODEL 

We consider the composite system of two coupled networks (I and II) with equal number 
N of oscillators for each. The equations of motion for the Kuramoto model in the system 
are given by 

N N 

% = ^-Jsin(0;-0; i )-J I ^a^sin(05-0j), $ = cof- Jsm{<f>f-<f>))- J n £ a* sinty?-^ 1 ), 

i i 

(1) 

where <f>j is the phase of the jth oscillator and is its intrinsic frequency in the 

network I (II), assumed as an independent quenched random Gaussian variable with the unit 
variance. Throughout the paper, we denote the ID regular network as the network I and the 
WS network as II, and use the average degree (k) = 4 for both. The adjacency matrix for 
the ID network has elements a* 7 = $j,i±i + $j,i±,2, while a5 is constructed following the small- 
world network generation method: Each link is visited and rewired at the probability p (see 
Ref. [13] for details). In Eq. fll]), J denotes the internetwork coupling strength between I and 
II, and Jmn stands for the intranetwork coupling for the network I(II). One may expect that 
when J becomes large enough the partial synchronization in the ID regular network could 
be induced by the established ordering in the WS network although no spontaneous ordering 
exists in the pure ID network. On the other hand, if J\ is not big enough, the internetwork 
coupling J to the nonsynchronized ID network could also make the WS network itself hard to 
synchronize. It is also possible that if J\ is infinite, yielding fully synchronized ID oscillators, 
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synchronization is induced in the WS network even at J\\ = 0. The interplay among these 
inter- and intra-network couplings can provide rich phenomena, the understanding of which 
composes the main motivation of the present study. 



III. ANALYTIC RESULTS 

A. Mean-field analysis for J = 

When the two networks are decoupled, i.e., when J = in Eq. 0]), the network II should 
show the MF synchronization transition as reported in previous studies |2J. In this Section, 
we briefly review the MF theory of the Kuramoto model for the single WS network 
The equations of motion for the network II without J is rewritten as 



JnAftfjSinW.' 1 -^), (2) 



where we have used Hje l9j = ct^e 1 1 with k^ 1 being the degree of the jth node 

in II. In the spirit of the MF approximation, we neglect the fluctuation and substitute Hj 
and 8j by global variables H and 9, respectively, which yields the self-consistent equation 
for the order parameter H: 



1 N ( to 11 \ 2 / \u u \ 



(3) 



where Q(x) is the Heaviside step function [Q(x) = 1(0) for x > (< 0)]. In thermodynamic 
limit of iV — > oo, we change the sum over oscillators to the sum over different degrees, which 
gives us 

H = -jz-^y^pikujknuiknHJn), (4) 
where p(kn) is the degree distribution function and 



u(x) 



/ duj n g{uii)^Jl - ulJx 2 (5) 

J — X 

with g{oJi\) = exp(— cjj I /2(Tj I )/-\/27ro"ii. Note that near the critical point where H becomes 
vanishingly small, Eq. ((5]) is expanded in the form u ~ |g(0)fcnif J\\ + f^g" (0)[knH Ju] 3 , 



which leads to 



Here, it is to be noted that since the network II is the WS network with the exponential 
degree distribution, both (kjj) and (k^) have well-defined finite values. It is then straightfor- 
ward to get the critical point J* 1 = ^ry^=— and the critical exponent f3 = 1/2. Moreover, 
introducing a sample-to-sample fluctuation [9| to the right-hand side of Eq. (jH]), we obtain 
the finite-size scaling (FSS) exponent v = 5/2. 



B. Mean-field analysis for J\ <C 1 

We next turn our attention to the coupled system (J / 0), and first consider the case 
of vanishingly small intranetwork coupling for I. As J\ — > 0, dynamics in the network I is 
simply governed by <fi\ = u>i — Jsin(0i — 0n) with the site index omitted for convenience. 
For (ui) = (uJii) = 0, running oscillators in I with <pi ^ may hinder their connected 
counterpart oscillators in II from entering into the global entrainment. It has also been 
found that contributions of detrained oscillators to the synchronization order parameter are 
negligible within the MF theory js, sj. Accordingly, one can make the plausible assumption 
that entrained oscillators (0i = 0) in I and their corresponding oscillators in II mainly 
contribute to the synchronization. We then write the equations of motion for the entrained 
oscillators in II as 

0n = w n +&i- JuhiH sin(0 n - 0), (7) 

where ui = Jsin(0i — 0n) with < J, and the MF approximation (Hj,6j) — > (H,6) has 
been made in the assumption that the internetwork coupling does not change the universality 
of II since there exists no ordering in I for J\ — > 0. Consequently, the self-consistent equation 
reads H = (iV'^n))" 1 £f kfyjl-ffe(l - |/,|), where f 3 = (uf + u})/kfHJ n and N> = 
N J^j duig(ui). In thermodynamic limit of N — >■ oo, we again meet the form of H = 
(V (hi)) J2k n p{hi)hiu(x) with 



/J rx—U>i j / _|_ \ 2 

duig(uji) J duj u g(ujii)J 1 - — ~J , (8) 

where x = kuHJu and g(ui) = J\T^ l e~ U) ^ 2(T ^ with the normalization constant J\f = 
J_j duie^^ / 2u i . Following the similar steps to those made for J = 0, we conclude that 
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the synchronization transition occurs at the critical value J* 1 = w ith a constant 



(9) 



We note from the expression that as J is increased, J* 1 also increases since A is decreased 
due to the fact that °/ P\ < of. Since a\ = crn, one obtains J 1 } {J) = tctt^ = a/2J"(0) 
with J — )• oo. Furthermore, it is expected that as J\ is increased (but still J\ <C 1) J* 1 should 
decrease from the following reasoning: When Jj 1, effective equations for oscillators in I 
can be written as (pi = oj[ — J sin(0i — n )- Here, u{ is a modified frequency whose variance 
a[ becomes smaller than the bare value due to the attractive force activated by the existence 
of the intranetwork coupling. We then conclude that a[ = cr[(Ji) is a decreasing function of 
J\. Substituting a\ by a[ in Eq. ([9]), one notes that r(Ji) = o^jio'i + ofj) increases with Jj, 
and thus A = (y^iwf^ 1 ff d dun exp [-w?/ (2ro"{ 2 )] with A/" = do^ exp [- Wl 2 / (2^ 2 )] 
becomes an increasing function (i.e., J* 1 becomes a decreasing function) with respect to J\. 
In summary of this subsection, when the ID network is within the weakly coupled regime 
(i.e., when J\ is sufficiently small), the synchronization of the WS network is enhanced (i.e., 
J* in reduced), as J\j J is increased. In words, the stronger internetwork coupling puts more 
burden for the WS network to achieve its synchrony, while the better synchrony in the ID 
network helps the WS network to be better synchronized via the internetwork coupling. 



C. Renormalization group approach for J\> J 

For a strong coupling regime of the ID regular network, our MF equations in Sec. IIIIBI 
for J\ — > need to be modified by employing the real-space renormalization-group (RG) 
formulation that has been developed in the ID systems [19[. In this RG approach applied 
for ID regular networks, strong bonds form clusters of entrained oscillators, while fast moving 
oscillators are decimated, interrupting the development of a giant synchronized cluster. In 
our system, it is expected that the strong coupling in the ID regular network should induce 
synchronized clusters not only in the ID regular network, but also in the WS network coupled 
to it via the internetwork coupling. 

For J\ > J, oscillators in the network I are governed mainly by intranetwork coupling 
rather than by internetwork coupling, which allows us to apply the RG approach to the ID 
network with respect to uj and J\. For > J\ > J, oj\ becomes relevant, whereas for 
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J\ > \u>i\ > J and J\ > J > \ui\, J\ becomes relevant. Applying the RG approach, fast 
moving oscillators having \ui\ > J\ are removed from the ID regular network together with 
their bonds, yielding fragmentation in I (our numerical RG calculations are summarized in 
Sec. UVp . On the other hand, since the remaining bonds are strong enough, oscillators left in 
a fragment form a synchronized cluster with the phase and the renormalized frequency 
ies 1 wnere m j = is a number of oscillators in the jth cluster Sj. 
The entrainment of oscillators into clusters in the network I tends to drive oscillators in 
network II to cluster themselves in the same corresponding clusters because of the ferro- 
magnetic intercouplings. For simplicity, we assume that the same number N c of clusters are 
generated both in I and II. Within the framework of the MF approximation, equations of 
motion are written as 



$i = n) - Jsin($5 - <S>J), $f = nf - Jsin($; 1 - $;.) - (k n )J n Hsm($f - 6), (10) 

where fi? = (1/m* 1 ) Yli^s 11 u l with being the number of oscillators in the jth cluster S 1 ] 1 
in the network II. A set of oscillators {</>* } connected to are merged into except for ones 
with | a;ii| > J since we only consider entrained clusters induced by the internetwork coupling. 
(&n) in Eq. (TTU|) also arises from (l/?7i*- ) ^2 ie gu k] ~ (kn). Now, it is straightforward to 
obtain J* 1 from Eq. ffTUj) via the same procedure as done for the case of J\ — 0, to yield 



where the approximation J^j dtti « d£li has been made from the assumption that 
o"i < J. Here, ofm) is the variance for fii(ii) within the Gaussian approximation, that is, 
g(Q m ) = e~ n i( II ) /2 ^(ii)/y / 27r^J". Note that a x ~ cr^/y/mi and a n ~ aJ/Vmi, where 
mi(n) is the average number of oscillators forming a cluster in I(II), given by m^n) = 



N f-}i(-j) dui(n)g(ui(ii)) / N c and cr* (n) is obtained from a^ 1T) = f^Lj) ^i(ii)^i(n)^(^i(ii)) 
since the frequency of oscillators which consist of clusters in I(II) must be less than J\ (J). 
For sufficiently large J and J\ yielding cr^ n ^ ~ 1 and rni(n) ~ N/N c , it is obtained that 
J* 1 ~ a/ N c /N , implying that the synchronization onset decreases with J\ since the number 
of clusters should be the decreasing function with large Jf, only at J\ — > oo, N c /N becomes 
close to 1/N, yielding J* 1 = in thermodynamic limit. 



*2 _ rMJ) a 2 
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We have above investigated the synchronization onset as a function of J and Ji via the MF 
approximation. It has been found that for weakly coupled ID oscillators, the internetwork 
coupling increases the synchronization onset, while increasing Jj enhances the synchroniz- 
ability. The role of the strong intranetwork coupling in the ID regular network has been 
revealed through the RG approach; clustering in the ID network makes the synchronization 
easier. 

For FSS exponent for J\ > 0, v = 5/2 still holds since the number of samples, i.e., N' or 
N c is proportional to the system size N linearly: The density N c /N of clusters only depends 
on Ji since oj\ is a random variable, yielding N c ~ N. 



IV. NUMERICAL RESULTS 



In order to validate our MF predictions, we perform numerical integrations of Eq. ([T]), 



13| at the rewiring rate p = 0.5. After 



with the network II constructed via the WS model 
achieving the steady state, we measure the time-averaged order parameter 

N 



A 



i(ii) 




E 



(12) 



where (...) stands for the average over different realizations of frequencies and networks. 

First, we compute order parameters when the ID network is in the weakly coupled regime 
(small Ji) at fixed Jn. In Fig. |2j we exhibit An and Aj as functions of J\ in (a) and (b), 
respectively. It is observed that as the internetwork coupling J is increased, the order 
parameters Ai and An show a decreasing tendency, except for very small values of J for Ai. 
The nonmonotonic change of Aj with respect to J is not surprising, since Aj ~ as J — > 
[see the curve for J = 0.1 in Fig. EJ^b)]. 

For Ji = 0, the synchronization onset J^ 1 is evaluated using the FSS of the form 



A n = N-WfdJn- J c U )iV 1/p ) 



(13) 



with p = 1/2 and v = 5/2, resulting in J c n (J) « 0.66, 0.68, 0.74, 0.8, 0.88, 0.89 and 
0.91 at J = 0, 0.5, 1.0, 1.5, 2.0, 3.0 and 4.0, respectively. In Fig. El it is observed that 
the ratio J ( [ I (J)/J^ I (0) obtained numerically follows the MF prediction in Sec. IHI B[ given 
by J,J I (J)/J ( [ I (0) = Jq 7 dxe~ x2 / 2 / dxe~ x2 with o\ = an = 1. It is also found that J* 1 
saturates at large J: J ( [ I (4.0)/J ( ! I (0) ~ 1.4, which is consistent with the MF value, J<F(J — > 
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FIG. 2: (Color online) Order parameters obtained through the numerical integrations of Eq. (fTJ): (a) 
for the WS network (An) and (b) for the ID regular network (At) versus the intranetwork coupling 
strength in the ID network Jj, at various values of internetwork coupling J (the intranetwork 
coupling Jjj in the WS network is set to 0.7 and the network size N = 12800 for both I and 
II). For sufficiently small values of Jj, the synchrony in II is suppressed (enhanced) as J(Ji) is 
increased, as shown Sec. IIIIBI via the MF approach. Inset in (a): Further increase of J\ enhances 
the synchronization in II. The internetwork coupling J plays a positive role for the synchronization 
in II, differently from the weak-coupling regime: As J is increased at a fixed Ji, An is increased, 
(b) shows that for small values of J\ the increase of internetwork coupling first enhances synchrony 
in I and then reduces it as J is increased further. 

oo)/J ( ! I (J = 0) = a/2. Also in the weakly-coupled regime of small Jj(> 0), it is again 
observed that J* 1 is increased as J is increased, implying that the stronger internetwork 
coupling worsens the synchronizability of II, in agreement with the finding in Sec. IIIIBI 
Another MF prediction in Sec. IIII Bl that the increase of J\ enhances the synchrony in II 
is clearly confirmed in Fig. |2](a): At fixed J, the order parameter for II is an increasing 
function of J\. For larger Jj_, the J-dependence of An shown in the inset for Fig. [2fa) can be 
interpreted as follows: More oscillators should be involved in synchronization with stronger 
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FIG. 3: (Color online) Ratio J c n (J)/J c n (0) as a function of J at Jj = 0. Points are from the 
numerical integrations of Eq. (pQ) combined with the finite-size scalings in Eq. (|13p , while the line 
denotes the MF prediction made in Sec. IIIIBi 

internetwork coupling in I, and thus the order parameter increases as J is increased, which 
may be implying that J* 1 is decreasing with J at large J\ in contrast to the case for the weak 
coupling regime of small J\. Consequently, we conclude that the role of the internetwork 
coupling J is reversed in the weak and the strong coupling regimes: For the weakly (strongly) 
coupled ID network, internetwork coupling worsens (enhances) the synchrony in the WS 
network. 

Very interestingly, a reentrant behavior of the order parameter as a function of Jj is 
observed at certain values of J. For example, at J = 0.5 and Jq = 0.55, nonmonotonic 
behaviors of An are displayed in Fig. H] (a): An(Ji) increases with J\ in accord with the 
prediction in Sec. IIII Bl begins to decreases at around J\ ~ 1.2, and finally increases again. 
The increase of An at large values of J\ is consistent with the calculation made in Sec. IIII CI 
In Fig. H] (a), also shown is that the nonmonotonic behavior does not change much with 
the system size. We observe that this reentrance is seen only in the limited range of the 
internetwork coupling, and disappears at lager J, as displayed in Fig. 0] (a) (see the upper 
most curve for J = 1.0). In Fig. H] (b) and (c), we exhibit scaling behaviors of An with N. 
What is found is that MF-like critical behavior, An ~ N~PI U with f3/i> = 1/5, is seen both 
at Ji ~ 1 and J\ ~ 4, whereas An decays more rapidly with iV at other values [e.g., J\ = 0.4 
and 2.0 in Fig. 11(b) and (c)]. We believe that this observation clearly indicates that as J\ 
is increased the curve of A n (Ji) touches phase boundary twice, first at J\ w 1, and later at 
Ji 4, implying the existence of a reentrant transition. 

To concrete our conclusion, we perform the FSS for An as a function of Jn with varying 
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FIG. 4: (Color online) (a) An versus J\ at J = 0.5 and Jn = 0.55 for various system sizes. 
Nonmonotonic change of An is observed at all system sizes, ^-dependence of An for (b) J\ =0.4, 
1.0 and (c) 2.0, and 3.8 at J = 0.5 and J n = 0.55. The upper (lower) line is for N' - 2 (A^ - 5 ). 
The MFdike scaling N~ ' 2 is found at J\ ~ 1 and at J\ ~ 4, indicating that the system undergoes 
MF transition twice as J\ is varied. In (a), for comparison, An for J = 1.0 and Jn = 0.5 (the curve 
at the top) is shown to be monotonic. 

J h as shown in Fig. El Since N^ p A n = f (0) at Jn = J c n from Eq. crossings shown in 
Fig. |5] clearly manifest the reentrance behavior of J* 1 : The crossing point decreases as Ji is 
increased from 1.0 to 1.3, and then increases as J\ is increased to 2.0. We report J* 1 obtained 
from FSS for J = 0.5 in Fig. [6j J l }{J\) first decreases from J\ = as predicted in Sec JIIIB[ 
and makes an up turn before eventually decreases again for larger Jt, consistent with Fig. HI 
Here, we emphasize that the reentrance behavior of the onset develops at J\ > J. 

The hump structure of the synchronization onset at J\ > J in Fig. [6] can be understood 
qualitatively by recalling the MF theory for J* 1 at Ji > J [see Eq. (jTTj) ]. given by a function 



of fluctuations of frequency in clusters, J^{J\) ~ a/ of + <x n . Rewriting Eq. ffTTj) as 

J c n (Ji) ~ J(N e /N) (A/f Vf 2 + A/^Vfj 2 ), 



with (Tj/jjn 



(14) 

•^1(11) I-i^-j)^m) u \ii)9{^i{U)) and TVi(ii) = ^(n^O^n))* one c& n 

find ingredients that determine the onset: For fixed J, A/'jjVj 2 is constant, while A/^Vj 2 is 
an increasing function of J\. It is expected that as J\ is increased, N c increases at small J\ 
since clusters are created rather than merged. For large Jt, on the other hand, N c should 
decrease as distinct clusters are merged. 
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FIG. 5: (Color online) Finite-size scaling [see Eq. (j!3j) ] for An at J = 0.5 and Ji = (a) 1.0, (b) 
1.3, and (c) 2.0. Crossing point first decreases and then increases, indicating that there exists a 
reentrant transition. 




FIG. 6: (Color online) Phase boundary at J = 0.5 separating the synchronous (the upper region) 
and the asynchronous (the lower region) phases. Points are obtained from the FSS of An [see 
Fig. [5]. As Ji is increased from below at J\\ ~ 0.56, the system starts from asynchronous phase, 
enters the synchronous phase, and then reenters back to the asynchronous phase as noted by the 
horizontal arrow. 
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FIG. 7: (Color online) Numerical RG calculation of N c /N as a function of J\ in the ID regular 
network with (k) = 4. One thousand different sets of intrinsic frequencies are used for the average. 

We finally perform numerical RG analysis for the ID regular network as described in 
Sec. IIII Ci Initially, the Gaussian frequency is distributed for a ID network of the size N 
up to 10 6 , and oscillators with \cjj\ > Ji are removed together with their bonds. Collecting 
fragments of the network, we calculate N c as a function of J\. From the numerical RG 
analysis, we find that N c /N as a function of Jj exhibits a well-defined peak near Jj f» 0.45 
as seen in Fig. [7J This allows us to expect that Jl l {Ji) has a peak at Jj* slightly larger 
than 0.45 since of in the right-hand side of Eq. (1141) is also an increasing function of J\. 
Since the expression for J* 1 in Eq. f fl~4"l) is valid only for J\ > J, for larger J than J x *, the 
reentrance behavior disappears as J is increased further, as seen in Fig. HI Although the 
values of J\ where the hump structure develops are different from those from the numerical 
investigation, the MF result made in Eq. ( fl~4l) gives us a qualitatively correct prediction. It 
indicates that when N c increases with increasing Ji, the synchronization in the WS network 
becomes worse since ID clusters with different frequencies would act as burdens, while the 
merging of clusters yields the better synchronization for sufficiently large J\. 

V. SUMMARY 

We have investigated synchronization phenomenon in the interdependent two coupled 
networks; one is the ID regular network and the other is the WS small-world network. Both 
the mean-field approximation and the numerical simulations have shown that the effect of 
the internetwork coupling is two folds: it suppresses the synchronization in the WS network 
when the ID network is in the weak-coupling regime, while it enhances the synchronization 
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in the strong-coupling regime. In comparison, the intranetwork coupling in the ID network 
has been shown to always play a positive role in the synchronizability of the WS network. 
In the intermediate range of the intranetwork coupling, the reentrant behavior has been 
found numerically and explained within the MF scheme combined with the numerical RG 
calculations. 

Acknowledgment 

This work was supported by NAP of Korea Research Council of Fundamental Science & 
Technology. B.J.K. and J.U. thank the members of Icelab at Umea University for hospitality 
during their visit, where this work was initiated. 



[1] A. T. Winfree, The Geometry of Biological Time (Springer- Verlag, Berlin, 1980). 

[2] A. S. Pikovsky, M. Rosenblum, and J. Kurths, Synchronization: A Universal Concept in 

Nonlinear Science (Cambridge University Press, Cambridge, 2001). 
[3] S. Boccaletti, V. Latora, Y. Moreno, M. Chavez, and D.-U. Hwang, Phys. Rep. 424, 175 

(2006). 

[4] Y. Kuramoto, in Proceedings of the International Symposium on Mathematical Problems 
in Theoretical Physics, edited by H. Araki (Springer, Berlin, 1975); Chemical Oscillations, 
Waves, and Turbulence (Springer, Berlin, 1984). 

[5] L. M. Pecora and T. L. Carroll, Phys. Rev. Lett. 80, 2109 (1998); M. Barahona and L. M. 
Pecora, ibid. 89, 054101 (2002); H. Hong, B. J. Kim, M. Y. Choi, and H. Park, Phys. Rev. E 
69, 067105 (2004); S. M. Park and B. J. Kim, ibid. 74, 026114 (2006); T. Nishikawa and A. 
E. Motter, ibid. 73, 065106(R) (2006); J. Um, B. J. Kim, and S.-I. Lee, J. Korean Phys. Soc. 
53, 491 (2008); S.-W. Son, B. J. Kim, H. Hong, and H. Jeong, Phys. Rev. Lett. 103, 228702 
(2009) ;R. T6njes, N. Masuda, and H. Kori, Chaos 20, 033108 (2010). 

[6] J. Gomez-Gardenes, Y. Moreno, and A. Arenas, Phys. Rev. Lett. 98, 034101 (2007); Phys. 
Rev. E 75, 066106 (2007). 

[7] H. Hong, M. Y. Choi, and B. J. Kim, Phys. Rev. E 65, 026139 (2002). 



15 



[8] J. G. Restrepo, E. Ott, and B. R. Hunt, Phys. Rev. E 71, 036151 (2005); Chaos 16, 015107 
(2005). 

[9] H. Hong, H. Park, and L.-H. Tang, Phys. Rev. E 76, 066104 (2007). 
[10] H. Hong, H. Park, and M. Y. Choi, Phys. Rev. E 70, 045204(R) (2004); 72,036217 (2005). 
[11] H. Hong, H. Chate, H. Park, and L.-H. Tang, Rhys. Rev. Lett. 99, 184101 (2007). 
[12] P. Erdos and A. Renyi, Publicationes Mathematicae Debrencen 6, 290 (1959). 
[13] D. J. Watts and S. H. Strogatz, Nature (London) 393, 440 (1998). 

[14] S. V. Buldyrev, R. Parshani, G. Paul, H. E. Stanley, and S. Havlin, Nature 464, 1025 (2010). 
[15] For more examples, S. Havlin, N. A. M. Araujo, S. V. Buldyrev, C. S. Dias, R. Parshani, G. 

Paul, and H. E. Stanley. larXiv:1012.0206V l. 
[16] R. Parshani, S. V. Buldyrev, and S. Havlin, Phys. Rev. Lett. 105, 048701 (2010); Proc. Natl. 

Aca. Sci. (USA) 108, 1007 (2011). 
[17] H. -H. Jo, S. K. Baek, and H. -T. Moon, Physica A 361, 534 (2006). 

[18] P. A. Robinson, C. J. Rennie, and D. L. Rowe, Phys. Rev. E 65, 041924 (2002); J. D. 

Victor, J. D. Drover, M. M. Conte, and N. D. Schiff, Proc. Natl. Aca. Sci. (USA) (in press) 

(doi:10.1073/pnas.l012168108). 
[19] O. Kogan, J. L. Rogers, M. C. Cross, and G. Refael, Phys. Rev. E 80, 036206 (2009); T. E. 

Lee, G. Rafael, M. C. Cross, O. Kogan, and J. L. Rogers, ibid. 80, 046210 (2009). 



16 



